Azbel-Hofstadter model on triangular lattice revisited 



O 

o 

(N 
5h . 

< 
(N 



Min-Fong YangQ 
Department of Physics, Tunghai University, Taichung, Taiwan 
(February 1, 2008) 

In the present paper, the mean of Lyapunov exponents for the Azbel-Hofstadter model on the 
triangular lattice is calculated. It is recently proposed that [Phys. Rev. Lett. 85, 4920 (2000)], for 
the case of the square lattice, this quantity can be related to the logarithm of the partition function 
of the two dimensional Ising model and has a connection to the asymptotic bandwidth. We find that 
the correspondence of this quantity to the logarithm of the partition function of the two dimensional 
Ising model is not complete for the triangular lattice. Moreover, the detailed connection between 
this quantity and the asymptotic bandwidth is not valid. Thus the conclusions for the mean of 
Lyapunov exponents suggested previously depend on the lattice geometry. 

PACS number(s): 73.43.-f, 71.20.-b, 71.10.Pm 
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The problem of two-dimensional (2D) Bloch electrons 
in a perpendiculpjpfliagnetic field has long been a fas- 
cinating subjecting Besides exhibiting extremely rich 
energy band structures, it can be related-to various phe- 
nomena such as the quantum Hall effejd.,0-3 the flux-state 
model for high-T c superconductivityJ13 and the mean- 
field transition temperature otsuperconducting networks 
or Josephson junction arrays£j In addition to the theo- 
retical work on this subject, due to the rapid develop- 
ment of nanofabrication techniques, the experimental in- 
dications |-efpthc peculiar band structure have also been 
reported.M 

The simplest model for this investigationjs-|Commonly 
referred to as the Azbel-Hofstadter model!!™ Although 
much work has been carried out, the analytic description 
of the spectrum remains open to a large extent. Recently, 
some progresses are made to relate this model to the in- 
tegrable systems, and thus deepen our understanding. 

It was noticed that the model can be related-to inte- 
grable systems with a quantum group algebra.Em H ence , 
explicit solutions to this problem can be constructedilMHl 
Recently, Krasovsky shows that a particular quantity, the 
mean of Lyapunov exponents, plays an important role in 
the study of the Azbel-Hofstadter problem and .gives a 
different type of relation to integrable systems £3 It is 
found that, for any value of magnetic flux $ = 2irP/Q 
through a plaquette (P and Q are mutually prime inte- 
gers), this quantity can be related to the logarithm of 
the partition function of the 2D Ising model. Moreover, 
the band edges of the Azbel-Hofstadter Hamiltonian cor- 
respond to the critical temperature of the Ising model. 
Besides the relation to the solvable 2D Ising model, it 
is proposed that the mean of Lyapunov exponents can 
also be connected to the asymptotic (large Q) bandwidth. 
These results provide a very exciting connection between 
two important models, and will create a perspective for 
even deeper understanding, if they are universal irrespec- 
tive of various kinds of lattice structures. For example, it 
is well-known that the concept of universality was discov- 
ered through the fact that the critical exponents of phase 
transitions in the Ising model are independent of lattice 



geometry. Thus it is interesting to investigate if these 
conclusions are still valid for other kinds of 2D lattice 
structures. 

The aforementioned theoretical and experimental stud- 
ies focused on the energy spectrum of a lattice with a 
square symmetry. Triggered by fascinating results for 
the square lattice, other work has considered various 
kinds of 2D lattices in order to investigate the effect 
of the lattice geometry on the energy spectrum. Be- 
sides the square lattice, the most-extensively studied case 
has been the triangular latticeE3~E3 The energy spec- 
trum of the triangular lattice with isotropic hopping is 
known to exhibit a recursive band structure as in the 
case rfrtfre square lattice.cia Moreover, it has been shown 
that£§EJ the same as it is for the square lattice, the 
sum of all bandwidths W ~ 32G/irQ as Q — > oo, where 
G = 0.91596559417721... is Catalan's constant. (Nota- 
tion A ~ B means, henceforth, that A/B tends to 1 in 
the limit.) However, in contrast to the case of the square 
lattice, it is shown that, for the case of the triangular 
lattice, not only the symmetry of the energy spectrum is 
lowered, but also a degeneraewnat-fche band center of the 
energy spectrum is removed. EaEZtEa 

Motivated by Krasovsky's workjll in the present pa- 
per, we calculate the mean of Lyapunov exponents for 
the case of the triangular lattice, and investigate if 
Krasovsky's conclusions will depend on lattice structures. 
For simplicity, we assume isotropic hopping and set the 
hopping amplitude to be unity. We find that the mean 
of Lyapunov exponents of the Azbel-Hofstadter model 
still has a similar form to the logarithm of the partition 
function of the 2D Ising model. However, in contrary 
to the case of the square lattice, we show that the band 
edges of the Azbel-Hofstadter model do not always cor- 
respond to the critical temperature of the Ising model in 
the present case. It means that the identification between 
these two models is not complete in the case of the tri- 
angular lattice. Later, another integral form of the mean 
of Lyapunov exponents is derived. From this integral 
representation, it is easily found that the mean of Lya- 
punov exponents at zero energy is positive, rather than 
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zero as in the case of the square lattice. From this ob- 
servation, we conclude that, while this quantity and the 
asymptotic bandwidth both scale as 1/Q, a detailed re- 
lation between these two quantities, which is conjectured 
by Krasovsky, fails in the case of the triangular lattice. 
Thus the aforementioned role played by the mean of Lya- 
punov exponents in the study of the Azbel-Hofstadter 
problem depends on the the lattice geometry. 

In the investigation of the Azbel-Hofstadter problem, 
there are two complementary, but mathematically equiv- 
alent, approaches: one can either consider the strong 
magnetic field limit with, the lattice potential treated 
as a small perturbation,H or study the influence of a 
weak magnetic field on the strong lattice potential in 
the tight-binding approximations Here we take the latter 



approach. For the tight-binding case, the model on the 
triangular lattice is topologically equivalent to that on 
the N x N square lattice with-tha iiext-nearest-neighbor 
hopping in only one direction£j~E3'EI When the magnetic 
flux through a plaquette is 2irP/Q, the quasi- momenta 
k° x = 2irm/N and k y = 2-nk/N arc defined in the mag- 
netic Brillouin zone: < fc° < 2-k/Q and < k v < 2tt 
(that is, m = 0, 1, . . . , N/Q - 1 and k = 0, 1, . . . , N - 1). 
The matrix of the Azbel-Hofstadter Hamiltonian H can 
be decomposed into a direct sum of the Q x Q matri- 
ces Hkm for each fc and m. In the case of the triangular 
lattice, the dependence of the spectral determinant (char- 
acteristic polynomial) of H km on k and m can be shown 
as (for example, see Eq. (4.1) in Ref. p7|) 
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where cr(e) = £® + . . . is a polynomial of Qth degree whose coefficients do not depend on k and m. Therefore, 
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In the limit of infinite lattice (N — » oo), one can replace 2-KmQ/N and 2nkQ/N by continuous parameters x and 
y, respectively. As discussed in Ref. [l3[ the mean of Lyapunov exponents 7_ for the case of Q + P = odd integers 
can be written as 

7-(<r)= lim -Lln\det(H-eI)\ 

N— >oo iv 
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dx 



dy\n |cr(e) — 2 [cosx + cosy + cos(a; + y) 



(3) 



For the case of Q + P = even integers, the mean of Lyapunov exponents 7 + can be defined in the same way. By using 
the periodicity of the cosine function, one can show that 7±(— |oj) = 7 T (|er|). Thus we need only to consider the case of 
7_ (i.e., the case of Q+P = odd integers). We note that the energy spectrum for Q+P — odd integers is determined by 
the secular equation, er(e) = 2 [cosx + cosy + cos(a; + y)]. Because — 3 < 2 [cosx + cosy + cos(.t + y)] < 6, following 
the same discussion in Ref. [l3[ one can show that the spectrum of H in the limit N — > oo consists of Q bands, which 
are just the image of the interval [—3, 6] under the inverse of the transform a = cr(e). In particular, the band edges 
e ei satisfy <j(e ei ) = —3 or 6. 

On the other hand, for the_2D ferromagnetic Ising model with isotropic interactions on a triangular lattice, the 
partition function Z satisfies:LJ 
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lim — ;r In Z — In 2 = — ^ / dx dy In I C 3 + S 3 — S [cos x + cos y + cos(x + y) 
n^oo N 2 8tt z Jo Jo 



(4) 



2 



Here C = cosh(2 J/T) and S = sinh(2J/T), where T is 
the temperature and J is the interaction constant. It is 
found that Eqs. (||) and (|^) have similar integral repre- 
sentations. Thus one may conclude that the connection 
between these two models still holds in the case of the 
triangular lattice. Nevertheless, in contrast to the case 
of square lattice, we find that the band edges of H do 
not always correspond to the critical temperature of the 
Ising model for the triangular lattice. We note that, be- 
cause (C 3 + S 3 )/S > 3, the integrand in Eq. @ can 
be divergent only at the critical temperature such that 
C 3 + S 3 = 35, and it diverges at the upper/lower limits 
of integration, x = 0, 2tt and y = 0, 2tt. For the Azbel- 
Hofstadter Hamiltonian, at the band edges e ei satisfy- 
ing cr(e e J = 6, the integrand in Eq. (||) also diverges at 
x = 0, 2n and y — 0, 2tt. Thus there is a correspondence 
between these band edges of H and the critical temper- 
ature of the Ising model. However, at the band edges e ei 
satisfying a(e ei ) = —3, the integrand in Eq. (||) diverges 
only at x = y = 2tt/3 or 47r/3. It means that the band 
edges e ei satisfying a(e ei ) = —3 have no correspondence 
to the critical temperature of the Ising model. There- 
fore, the identification between these two models on the 
triangular lattice is not complete. 

It is noted by Krasovsky that ,113 the mean of Lyapunov 
exponents is proportional to 1/Q (see also Eq. (||) in the 
present paper). Since the total bandwidth is shown to 
be of order 1/Q as Q — > ooJ3 it is reasonable to think of 
a connection between these two quantities. For the case 
of the square lattice, Krasovsky shows that the mean of 
Lyapunov exponents at the band edges is eight times less. 
than the estimated asymptotic (large Q) bandwidth Otl 
By defining W(x) to be the total length of the image of 
[0, x] under the inverse of the transform a = a{e), he fur- 
ther conjectures that W(a) ~ 47(er) for any coprime P 
and Q as Q — > oo. It is interesting to investigate if this 
conjecture is still valid for the triangular lattice. 

For further discussions, it is more convenient to recast 
Eq. (||) into another from. By combining the terms cos y 
and cqs£.t + y), and then taking the integrals over y, we 
obtain^ 
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Hence for —2 < a < 6, one has 
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Here v = ^a + S, and K(k) = J*' (1 - k 2 sin 2 x^^dx 
is the complete elliptic integral of the first kind. The first 
equation in Eq. (Q) is derived by employing the standard 
result 
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(C 2 - i)- 1/2 , C> i ; 
-(C 2 -!)- 1 / 2 , C<-i; (8) 
o, ICI < i, 



and theJast equation in Eq. (Q) is obtained by using the 
formulal^J 
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Similarly, for — 3<cr<— 2, we obtain 
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da 
1 



dx 



27T<9 Jcos-^l-VF+E+f) ^/{cOS ■! 
^ r cos _1 (l+ A /?+3-l-f ) 



<T\2 
2 ^ 



4cos 2 (f) 



2ttq y 



A' 



(l-»/)V(l- 1 /)(3+^) 

(1+^)^(1+7^(3^7) 



^(cosa;- f) 2 -4cos 2 (f) 
(11) 



TTQ (1 + V )y/{\ + V ){i- V ) 

From these formulas, the mean of Lyapunov exponents 
at the band edges can be numerically calculated with 
the results: 7-(<r = 6) = 5.075/7rQ and j-(a = —3) = 
2.030/ttQ. Thus 4 7 _(cr = 6)+4 7 _(cr = -3) ^ 28.42/ttQ, 
which is close to the estimated value of the total width 
Wr=*rW(-3) + W(6) ~ 32G/ttQ = 29.31/ttQ as Q -> 
This seems to support the validity of Krasovsky's 
conclusion in the case of the triangular lattice. Never- 
theless, because d r y-(a)/da > for — 2 < a < 6, we find 
from Eq. @ that 7_(<r = 0) is positive rather than zero. 
(One can show numerically that 7_(a = 0) = 1.347/7rQ.) 
On the other hand, it is obvious that VF(0), which is the 
total length of the image of a zero interval under the in- 
verse of the transform a — a(e), must be zero. Therefore, 
the detailed relation, W(a) ~ 47(cr), which is conjectured 
in Ref. [l3[ fails in the present case. 

In summary, the integral expressions of the mean of 
Lyapunov exponents for the case of the triangular lattice 
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is obtained in the present work. From these expressions, 
we show that it can not be properly related to the log- 
arithm of the partition function of the 2D Ising model. 
Moreover, the conjectured relation W(a) ~ 47(a) can 
not be true for the general types of lattice structures. 
Therefore, the mean of Lyapunov exponents plays a less 
important role in the study of the Azbel-Hofstadter prob- 
lem on the triangular lattice. 
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